Short GRB and binary black hole standard sirens as a probe of dark energy 
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Observations of the gravitational radiation from well-localized, inspiraling compact object binaries 
can measure absolute source distances with high accuracy. When coupled with an independent 
determination of redshift through an electromagnetic counterpart, these standard sirens can provide 
an excellent probe of the expansion history of the Universe and the dark energy. Short 7-ray 
bursts, if produced by merging neutron star binaries, would be standard sirens with known redshifts 
detectable by ground-based GW networks such as LIGO-II, Virgo, and AIGO. Depending upon the 
collimation of these GRBs, a single year of observation of their gravitational waves can measure the 
Hubble constant h to ~ 2%. When combined with measurement of the absolute distance to the last 
scattering surface of the cosmic microwave background, this determines the dark energy equation 
of state parameter w to ~ 9%. Similarly, supermassive binary black hole inspirals will be standard 
sirens detectable by LISA. Depending upon the precise redshift distribution, ~ 100 sources could 
measure w at the ~ 4% level. 



I. INTRODUCTION 

With the advent of the Laser Interferometer 
Gravitational- Wave Observatory (LIGO), we are on the 
verge of an era of gravitational- wave astronomy 0, 0- 
Among the most interesting expected sources for GW 
observatories are the inspirals and mergers of compact- 
object binaries. LIGO-II, a planned upgrade with ten- 
fold increase in sensitivity, can detect the inspirals and 
coalescence of stellar-mass binaries within several hun- 
dred megaparsecs, while the Laser Interferometer Space 
Antenna (LISA) can study supermassive binary BHs 
(M ~ 10 4 - lO 7 Af ) throughout the universe (z < 10). 

The idea of using GW measurements of coalescing bi- 
naries to make cosmologically interesting measurements 
has a long history. As originally pointed out by Schutz 
0, observation of the gravitational radiation from an 
inspiraling binary provides a self-calibrated absolute dis- 
tance determination to the source. Chernoff and Finn 
and Finn p| took advantage of this property to show 
how, by observing many inspiral sources, one can con- 
struct the distribution of observed binary mass and GW 
signal strength, and thereby statistically constrain the 
values of cosmological parameters. More recently, Holz 
and Hughes [|| have shown that LISA observations of 
well-localized supermassive binary black hole (SMBBH) 
inspirals allow cosmological distance determination with 
unprecedented accuracy, with typical errors < 1%. These 
GW "standard sirens" can precisely map out the expan- 
sion history of the Universe, offering a powerful probe of 
the dark energy. 

The utility of standard sirens for constraining dark en- 
ergy is quite similar to that of standard candles, like 
Type-la supernovae. One advantage of GW standard 
sirens is that the underlying physics is well-understood. 
The radiation emitted during the inspiral phase (as op- 
posed to the merger phase) is well described using the 



post-Newtonian expansion of general relativity for the 
BH binary Q. Hence some unknown systematic evolu- 
tion of the standard sirens over time, mimicking a dif- 
ferent cosmology, should not be of concern. Another ad- 
vantage is that GW observatories directly measure abso- 
lutely calibrated source distances, whereas Type-la su- 
pernova standard candles provide only relatively cali- 
brated distances. 

A major drawback of GW standard sirens is that, al- 
though the gravitational waveforms measure distance di- 
rectly, they contain no redshift information. To be use- 
ful as a standard candle, an independent measure of the 
redshift to the source is crucial. This can be determined 
through observation of an electromagnetic counterpart, 
such as the host galaxy of the source. Unfortunately, as 
GW observatories are essentially all-sky, they generally 
provide poor source localization, and the host galaxy is 
not always unambiguously identifiable |8( ■ In cases where 
source redshifts cannot be determined, the distribution of 
unlocalized events can be used to place statistical bounds 
on cosmological parameters However, in this paper 
we will focus upon GW sirens whose redshifts may be 
measured, as they can provide very tight constraints on 
cosmology. 

Because standard siren distances are absolutely cali- 
brated, even sources at low redshift (e.g. z < 0.2) can 
constrain dark energy. This may seem surprising, since 
at low redshifts the distance-redshift relation is well- 
described by a linear Hubble relation D = cz/Hq, in- 
dependent of dark energy parameters. As emphasized 
by Hu and Jain [j| and Hu however, absolute dis- 
tances to sources at low redshift tightly constrain dark 
energy, when combined with a determination of the abso- 
lute distance to the last-scattering surface of the cosmic 
microwave background. To understand this, note that 
cosmological distances are given by a redshift integral of 
the Hubble parameter, which in turn depends on the sum 
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of energy densities at each redshift: 



II. DISTANCE DETERMINATION FOR 
INSPIRALING BINARIES 
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Here fl m + Cl de +fl K = 1, H = 100 h km/s/Mpc is the 
Hubble constant today, and we have assumed a constant 
equation of state parameter w. If we assume a flat uni- 
verse (0^- = 0), then f^c = 1 — ^m, and the only param- 
eters describing the global expansion are h, Q, m , and w. 
Observations of the primary anisotropics in the cosmic 
microwave background (CMB) provide two constraints 
on these three parameters. First, the heights of the 
acoustic peaks determine the matter density (in g/cm 3 ), 
which fixes Q. m h 2 . Second, the angular scale of the peaks 
(their location in /-space) precisely measures the angular 
diameter distance to the CMB last-scattering surface, in 
Mpc. Absolute distances to low-redshift sources measure 
the Hubble constant h, which then allows all three pa- 
rameters to be determined Clearly the constraints 
we present would be substantially degraded if the curva- 
ture were not fixed to zero; see [ll|,[l2| for the prospects 
for precise constraints on curvature. 

In addition to the low redshift standard sirens, those 
at higher redshifts also help constrain dark energy, in the 
same manner as high-redshift standard candles. Holz and 
Hughes discuss how LISA observations of SMBBH in- 
spirals can help constrain cosmology. For a dark energy 
model which is not dramatically different from a cosmo- 
logical constant A, the interesting redshift range is when 
the dark energy density is significant (z < 1), although 
note that gravitational lensing degrades the constraints 
from the highest redshift standard sirens (or candles) 

As mentioned above, the gravitational waves from 
standard sirens measure source distances, but do not 
measure source redshifts. Some sort of electromagnetic 
counterpart associated with the merger event will gener- 
ally be required to use GW sources to determine cosmol- 
ogy. One potential class of GW sources guaranteed to 
have electromagnetic counterparts are short 7-ray bursts 
(GRBs). These sources are thought to arise in the merg- 
ers of neutron star (NS) binaries, and hence should be 
strong GW emitters in the frequency band accessible 
to ground-based observatories. The GRB counterpart 
to these GW source provides a precise sky localization, 
which is useful both for determining the redshift to the 
source, and for significantly improving the GW determi- 
nation of absolute distance. As we discuss below, short 
GRBs occur at a rate large enough for them to provide 
interesting cosmological constraints. 



In this section we briefly review how distances to in- 
spiraling binaries may be determined; see Ref. |14| for 
more detail. An inspiraling binary at direction ft on the 
sky, with orbital angular momentum axis X, generates 
gravitational waves with strain tensor 



h(t) = h + (t)e + + h x (t)e > 
where the basis tensors arc 
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Our convention is that n points towards the source, hence 
the waves propagate in the direction — n. We can express 
the amplitudes of the two polarizations h+(t) and h x (t) 
in the frequency domain as 

h + (f) = (l+v 2 )h (f), h x (f) = -2ivh (f) (7) 

where v = h ■ L is the cosine of the inclination angle of 
the binary, and 
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In this expression, D is the luminosity distance to the 
source, and M = (1 + z)[mim2] 3 ^ 5 / (mj + m 2 ) 1 ^ 5 is the 
redshifted chirp mass of the binary. The phase ^ is given 
by 
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where t c is the time at coalescence, and <p c is the orbital 
phase at coalescence. 

These expressions describe a binary's waves only in the 
Newtonian, quadrupole approximation — treating the bi- 
nary's kinematics as due to Newtonian gravity and using 
the quadrupole formula to estimate its GW emission. Be- 
cause the phase parameters are essentially uncorrelated 
from the amplitude parameters, this approximation is 
good enough to estimate the expected signal-to-noise ra- 
tio from a source, and provides a good estimate of the dis- 
tance measurement accuracy, but is not accurate enough 
to reliably model the detailed GW waveform 01 • Higher 
order post-Newtonian templates (see Ref. 7) for detailed 
discussion) should be sufficiently accurate, and are used 
for the actual data analysis. 
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Given h(t), the measured strain is given by 

h M (t) = h ab (t)d ab , (10) 

where the detector response tensor for an interferometer 
with arms I and rh is d = (l®l — rh®rh)/2. In the nota- 
tion of Ref. 01 1 a detector at colatitude and longitude 
(j> with orientation a has response tensor 

d = cos(2a)[e e ® + eg]/2 

- sin(2a)[e 9 (g) e e - (g) e^]/2. (11) 

To recap, the source parameters determining the mea- 
sured signal are distance D, chirp mass A4, coalescence 
time t c , coalescence phase <p c , source direction n, and 
orbital axis L. These are the 8 parameters to be deter- 
mined from the data timestream /im(0- If t ne detector 
has strain noise with spectral density Sh{f), then the in- 
cident strain is measured with signal-to-noise ratio SNR 
(assuming Wiener filtering): 



SNR 



\h M (m 

ShU) 



df. 



(12) 



The complicated angular dependence is hidden within 
the measured strain Km- This dependence can be made 
more explicit by rewriting the above equation as p| 
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where A = ^f^jm tt" 2 / 3 (GM/c 3 ) 5/6 c, F+ = e+d' 
F x = e*d ab , and 
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Here /i ow — 10 Hz is the frequency below which the detec- 
tors' sensitivities are badly degraded by ground motions. 
In the optimal case, the binary is face-on (v = 1) and 
directly overhead, so that Ff + F 2 =1. This gives 
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If instead we average over all sky positions and binary 
orientations, we find 
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where we have made use of {FT) = (F 2 ) = 1/5 and 
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note that face-on sources, whe n av eraged over all sky 
positions, have SNR a factor y / 5/4 ~ 1.12 larger than 
SNR ave . 

We can estimate how well the parameters p are mea- 
sured using the Fisher matrix 
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where di = d/dpi, and * denotes complex conjugation. 
Approximating the likelihood as 
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then the error on parameter pi is given by U (F" 1 )^. 
Prior constraints, or the constraints from multiple de- 
tectors are implemented by multiplying the respective 
likelihoods, which in this approximation reduces to sum- 
ming the respective Fisher matrices. In our calculations, 
we compute the partial derivatives numerically by finite 
differencing. 

In practice, the 'phase' parameters M., t c , and <p c are 
determined with exquisite precision. The 'amplitude' pa- 
rameters Z), i, and ft are determined less well, in large 
part due to parameter degeneracies. By using multiple 
detectors many of these degeneracies can be broken. For 
example, timing information from a network of detec- 
tors helps determine the source direction n. Similarly, if 
the detectors have different response tensors d, then the 
polarization of the GW signal may be measured, which 
constrains the orbital axis L (c.f. Eq. (JJJ)). 



III. GRBS OBSERVED BY GW NETWORKS 

Short GRBs are an extremely promising source of grav- 
itational waves 15]. These sources have been of great 
interest recently, due to the prompt localization of the 
events by the Swift 1 [H and HETE-2 2 0| satel- 
lites, allowing their detection in X-ray, optical, and radio 
frequencies. Particularly exciting has been the identifica- 
tion of several galaxies hosting short bursts 0, 0, ^| . 
While the nature of short GRBs is not yet known, a lead- 
ing candidate is the merger of neutron star binaries [2pJ , 
although other models have been proposed as well |21j . 
The detection or non-detection of GRBs in gravitational 
waves would of course be extremely useful |22| , for exam- 
ple confirming or refuting the NS-NS merger scenario, or 
determining the extent of collimation of the 7-ray emis- 
sion H3. 

Additionally, as mentioned above, short GRBs can also 
be very useful for determining the background cosmology 



Note that the SNR in the optimal geometry is a factor 
5/2 times larger than that for the average geometry. Also 



1 http:/ /swift. gsfc.nasa.gov/docs/swift/swiftsc. html 

2 http://space.mit.edu/HETE/Welcome.html 
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by acting as GW standard sirens. One immediate advan- 
tage offered by GRBs is that their bright electromagnetic 
emission allows a precise localization of the source on the 
sky, pinpointing the source direction h and lifting some of 
the degeneracies which limit distance determination. The 
extent of collimation in short GRBs is not well known, 
though indications of beaming are claimed in at least 
one short GRB so far . The theoretical expectation is 
that emission should be beamed preferentially along the 
orbital angular momentum axis where baryon loading is 
minimized. If this is the case, then we might expect short 
GRBs to be nearly face-on, v = n- L » 1. As can be seen 
from Eq. (JJJ this maximizes the amplitudes of both po- 
larizations and hence maximizes the SNR of the detection 
for a given source direction fi. In what follows, we will 
compute distance errors for two cases: (1) isotropic distri- 
bution of i, and (2) collimation, assuming an inclination 
probability distribution dP/dv oc exp(— (1 — v) 2 /2a 2 ) for 
a v = 0.05, corresponding to a roughly 20° jet angle. 

The expected chirp mass for GRBs, M. ~ 1.2M©, 
places them favorably in the frequency band accessible 
to ground-based GW observatories. Several such obser- 
vatories are now operating or are planned for construc- 
tion in the near future. LIGO is already operational, 
and its sensitivity should increase by an order of magni- 
tude in a planned upgrade (LIGO-II) |24| . A detector of 
similar scale, Virgo 25], is under construction in Italy, 
and there are plans for a similar detector, AIGO j2(J , in 
Australia. The locations and orientations of these obser- 
vatories are listed in Table |U The two LIGO detectors 
are oriented to have very similar response tensors, and 
therefore have limited ability to independently measure 
polarization (and hence inclination). Determining L will 
thus require the combination of LIGO with other obser- 
vatories. 

Henceforth, we assume that all four detectors will ob- 
serve GRB events; in subsequent work we intend to inves- 
tigate how the distance errors degrade if one or more el- 
ements of this network are removed. Preliminary results 
indicate that reducing the size of the detector network 
does not substantially degrade our ability to determine 
distance (aside from the loss in total SNR) assuming that 
we can set a prior on the beaming factor (and hence on 
the inclination angle). If we cannot set such a prior, then 
losing sites in this network badly degrades our ability to 
determine distance to these sources. We emphasize this 
point to highlight the importance of modeling bursts, and 
the importance of having widely separated GW detectors 
around the globe. 

Figure pjplots the noise spectral density forecasted for 
LIGO-II [24J. Projected noise curves for the advanced 
detector configurations are not yet available for Virgo or 
AIGO, so for simplicity we use the LIGO-II curve for 
all the observatories in the network. For comparison, 
we also show the sensitivity for the currently operating 
LIGO observatories. 

With the response tensors for the elements in our 
network, and their noise spectra, we can now compute 
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FIG. 1: Noise curve for the LIGO detectors, for initial (dot- 
ted) and advanced (solid) sensitivity. 



the Fisher matrices and parameter errors for GRBs as 
a function of distance and location on the sky. For 
convenience, when computing the Fisher matrix we re- 
place the parameter pair {D,v — h ■ L} with the pair 
{(1 + v) 2 /D, (1 — v) 2 /D}, to avoid singularities in the 
limit v — > 1 when v and D become degenerate 0]. An- 
other difficulty that arises in the face-on limit is that 
the position angle of L, denoted if) by Ref. 0], becomes 
meaningless as v — ► 1. Including it as a parameter would 
cause the Fisher matrix to become singular in the face-on 
limit; we circumvent this difficulty using singular value 
decomposition to invert the Fisher matrix, zeroing any 
eigenvalue whose magnitude is 10~ 10 times that of the 
largest eigenvalue. 

Because the antenna response of each detector varies 
strongly with source direction ft, the parameter errors 
at any given distance D also depend strongly on n. We 
are interested only in average errors as a function of D; 
hence, for each D we average over 100 different orienta- 
tions L and n. For example, Fig. |2 shows the expected 
constraints for sources at distance D — 250 Mpc. Note 
that the errors significantly improve if it is assumed that 
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LIGO (Hanford) 


43.54 


-119.4 


171 


LIGO (Livingston) 


59.44 


-90.77 


243 


Virgo 


46.37 


10.5 


115.6 


AIGO 


121.4 


115.7 
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TABLE I: Coordinates of GW observatories, in the notation 
of Ref. |14| . All values are in degrees. 



5 



0.030 




0.000 



200 300 
D [Mpcl 



500 



FIG. 2: Distribution of measured distances for a source at 
D = 250 Mpc, averaged over 100 source directions n and ori- 
entations L. The solid curve shows constraints for randomly 
oriented sources, while the dashed curve shows constraints for 
collimated sources with a v — 0.05. 



sources are beamed towards us. 

Given the likelihood distribution dP/dD, we define 
the distance error as erf) = (D 2 ) — (D) 2 , where aver- 
ages are with respect to dP/dD. Figure plots <jd as 
a function of D. Our results appear roughly consistent 
with <jd/D oc D oc 1/SNR. Our best-fit linear scal- 
ing for unbeamed GRBs is <jd/D = D/(l.7 Gpc), and 
ctd/D = D/(AA Gpc) for collimation a v = 0.05. Hence- 
forth, we will assume these scalings when estimating cos- 
mological constraints from GW network observations of 
short GRBs. 



IV. COSMOLOGICAL CONSTRAINTS FROM 
STANDARD SIRENS 



As discussed in §1, a measurement of the Hubble con- 
stant h using GRBs, when combined with CMB con- 
straints, also constrains dark energy parameters. We use 
two measurements from the CMB: determination of the 
angular scale of the acoustic peaks I a, and determination 
of the matter density Sl m h 2 from the peak heights. Cur- 
rently, the WMAP satellite has measured I a = 300 ± 3 
and fl m h 2 = 0.14 ± 0.02 [27|. The error on Q m h 2 will 
soon decrease by a factor ~ y3 with the 3rd year re- 
lease of WMAP data. The Planck satellite is expected 
to measure fi m /i 2 to a fractional error of ~ 1%. 

The acoustic scale Ia = irD+fs*, where D+ is the dis- 
tance to the last-scattering surface at z — 1089, and s* 
is the sound horizon at decoupling, approximately given 
by = 144.4 Mpc {Vl m h 2 /0.14)" - 252 0. Given the 
dependence of these observables on the cosmological pa- 
rameters p — {h, Cl m , u>}, we can then estimate parame- 
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FIG. 3: Fractional distance errors as a function of source 
distance D. The + symbols are for unbeamed GRBs, while 
circles are for cr„ = 0.05. The two lines show the best-fit 
linear relations; note that there may be departures from linear 
scaling at the highest distances. 



ter errors using the Fisher matrix: 



dN diD L (z)djD L {z) 
dz a D {z) 2 + (a **2tV 



dz, (21) 



dz 



where redshift errors <j z are caused by peculiar velocities 3 
with assumed rms of 300 km/s. The luminosity distance 
Dl{z) = (1 + z)D(z), and its error an includes both GW 
errors, as computed in the previous section, and gravita- 
tional lensing errors 28] , computed using an approximate 
nonlinear power spectrum |29| . 

For the source redshift distribution dN/dz, we assume 
that short GRBs occur at a constant comoving rate of 
10 Gpc -3 yr _1 30]. We found in the previous section 
that the SNR in distance determination per source scales 
roughly like 1/D. Since the number of sources scales with 
volume oc D 3 , we expect the SNR on the Hubble constant 

1/2 

h to scale like -Dmax, where -D max is the maximum dis- 
tance to which GRBs may be detected as gravitational 
wave sources. 

The standard threshold used in the GW literature for 
detection has been SNR > 8.5 The reason for 

this high threshold is that sources are detected by cor- 
relating the data timestream with large numbers (e.g. 



3 It may be preferable to measure redshifts of the host galaxies 
rather than the GRBs themselves, whose progenitors may suffer 
kicks which will add in quadrature to the redshift noise from 
peculiar velocities. 
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FIG. 4: Errors on h and w as a function of GRB exposure, as- 
suming PiancJc-quality errors from CMB. Solid curves are for 
ah, the error on the Hubble constant, while dashed curves cor- 
respond to a w , for the dark energy equation of state param- 
eter. The lower curves are for beamed GRBs with a v = 0.05 
while the upper curves are for unbeamed GRBs. 



10 15 ) of templates corresponding to different parameter 
values, and therefore the detection threshold must be set 
high to avoid excessive numbers of false detections. Such 
large numbers of templates are required in order to fully 
explore parameter space. For GRB sources, however, the 
parameter space to be searched is considerably reduced: 
the 7-ray burst itself determines the source direction h 
and time t c . Depending upon one's confidence in theoret- 
ical models for GRBs, the chirp mass Ai and orientation 
L may also constrained. Because many fewer templates 
need to be run for GRB sources, we should set the detec- 
tion threshold correspondingly lower. We conservatively 
estimate that knowledge of the GRB time reduces the 
number of required templates by a factor ~ 10 5 , corre- 
sponding to a reduced threshold SNR > 7. Note that 
this is the total SNR; since we have assumed a network 
of four detectors with identical noise, this translates into 
a threshold SNR > 3.5 per detector. From this, we can 
determine the maximum distance to which sources may 
be detected using Eq. (jTo)l . For chirp mass M = 1.2M , 
we have A = 4.7 x 10 _6 s 5 / 6 , and for our assumed noise 
spectral density (Fig.Q), I 7 = 8.33 x 10 44 Hz~ 1/3 . There- 
fore the maximum distance for which SNR avc > 3.5 is 

Bmax = 600 M P C. 

Assuming default cosmological parameters h = 0.72, 
Q m = 0.27, and w = — 1, the resulting parameter er- 
rors computed from Eq. (|21|l are shown in Figure as a 
function of the time and sky area over which GRBs are 
observed. While errors on the Hubble constant scale like 

— 1/2 

cr/j cx -Nqrb ' t ne errors on w scale this way only in the 
limit of small numbers of sources. Quite rapidly, the lim- 
iting error on w becomes the uncertainty in CMB (in this 



figure, fractional errors of 1% on f2 m /i 2 were assumed). 
Unless CMB errors can be significantly improved, it will 
be difficult for low-redshift GW sources to constrain w 
to better than the ~ 10% level. 

Higher redshift standard sirens would probe departures 
of the cosmic expansion from linear Hubble scaling, and 
thereby directly constrain parameters like Q m and w. 
Unfortunately, stellar-mass inspirals at high redshift are 
not sufficiently luminous to be detected by any exist- 
ing or planned GW observatory. Inspirals involving su- 
permassive black hole binaries, however, are sufficiently 
luminous in gravitational waves to be detected at cos- 
mological distances. As discussed by Holz and Hughes 
, LISA observations of SMBBH inspirals can in princi- 
ple measure distances to better than 1% accuracy. This 
precision is degraded, however, by gravitational lensing 
caused by density fluctuations from large-scale structure 
along the line of sight to the source. Another difficulty 
in using LISA observations is that, unlike in the case of 
short GRBs, for SMBBHs there are no guaranteed elec- 
tromagnetic counterparts. However, it has been argued 
that many SMBBH mergers will be followed by bright 
quasar-like activity [32 > or possibly preceded by optical 
emission j33|, which will localize the GW source on the 
sky and provide a source redshift. 

Due to lensing errors, small numbers of LISA sources 
will generally be unable to constrain dark energy param- 
eters significantly [(| . The effects of lensing diminish sig- 
nificantly at lower redshifts, so a single SMBBH inspiral 
at z < 0.5 observed by LISA could measure the Hubble 
constant to < 1% and w to < 10%. Although such a 
source is unlikely, the low redshift regime should already 
be well-determined by ground-based GW observations of 
short GRBs. On the other hand, if large numbers of 
SMBBH mergers occur during LISA's lifetime, then LISA 
should provide quite interesting constraints on dark en- 
ergy, despite the lensing noise. To illustrate this, Figure[5] 
plots expected constraints in the Cl m vs. h plane for a 
sample of 100 SMBBH inspirals observed by LISA, dis- 
tributed in redshift assuming a constant comoving den- 
sity between < z < 2, combined with constraints from 
Planck-quality CMB data. The 1 — a errors on w are 
o~ w = 0.04; these are competitive with ambitious Type- 
la supernova surveys like JDEM. Note that these errors 
improve considerably if the main source of noise, grav- 
itational lensing, can be cleaned out by reconstructing 
the lensing mass distribution using other probes. Dalai 
et al. [25 argue that cosmic shear measured from optical 
surveys would not allow mass reconstruction with suffi- 
cient angular resolution. Cosmic magnification measured 
in the radio could conceivably offer an alternative route 
(e.g. HI). 

Our discussion has focused on gravitational lensing 
only as a source of noise, but in principle there is cos- 
mological information which can be extracted from the 
lensing fluctuations themselves. With large numbers of 
sources, LISA observations of cosmic magnification can 
provide constraints complementary to other probes. We 
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FIG. 5: LISA constraints on dark energy. The solid contours 
show the 68% and 95% confidence regions expected for a sam- 
ple of 100 SMBBH sources observed by LISA, distributed with 
constant comoving density between < z < 2. A Planck 
prior also has been used on fi m /i 2 and I a, as discussed in the 
text. The dotted contours correspond to a sample of 3000 SNe 
with intrinsic luminosity scatter of 10%, with redshift distri- 
bution oc exp(-(z - 0.5) 2 ) over 0.02 < z < 2. The dashed 
(dark shaded) contour shows the 68% confidence region for 
the combined constraints GW+SNe+CMB. 



would not expect GW standard sirens to usefully probe 
the power spectra of matter fluctuations or galaxy-mass 
correlations 9| at any scale, compared to other means 
like cosmic shear or Type-la supernovae, based on their 
noise power spectra : 



a powerful probe of cosmology. In particular, short 7-ray 
bursts appear quite promising as potential GW standard 
sirens. The presently observed rate of short GRBs is 
sufficiently high that within a few years of observation 
by the next generation of ground-based GW observato- 
ries (e.g. LIGO-II, Virgo and AIGO), strong constraints 
on dark energy parameters may be derived (cr w < 0.1). 
These inspiraling NS-binaries should be clean sources of 
gravitational waves; possible sources of contamination, 
such as tidal effects, magnetic torques, or gasdynamical 
torques from circumbinary gas, should all be negligible 
during the crucial inspiral phase (where v/c < 0.3). We 
emphasize again that the best information about distance 
measurements comes from combining multiple GW data 
from instruments that are widely separated. Good infor- 
mation about the collimation of the gamma rays and thus 
on the likely inclination of the binary progenitor will also 
improve the utility of these standard sirens. Given the 
great cosmological potential of GW observations of short 
GRBs, there is strong incentive to extend the lifetime of 
GRB satellites such as Swift or HETE-2 to overlap with 
next-generation gravitational wave observatories. 

The inspirals of SMBBH binaries observed by LISA 
can also provide interesting constraints on dark energy, 
if the rate of such mergers is high enough to average away 
noise caused by gravitational lensing. At present, the to- 
tal rate and redshift distribution of SMBBH mergers are 
not well understood, with estimates ranging from a few 
(or zero) per year, up to hundreds per year, depending 
upon assumptions (3(| 113, IH, ■ If the rates are at the 
high end of these estimates, with a significant fraction at 
redshifts z < 2, then w may be constrained at the few 
percent level. 



Tgal o 



2 

GW 



n ga i 



new 



(22) 



where galaxies have shape noise 7 ga i ~ 0.4 and num- 
ber density n ga i m 50/arcmin 2 , supernovae have lumi- 
nosity dispersion ctsn ~ 0.1 and number density nsN ~ 
4000/(20 deg) 2 as observed by SNAP, and GW standard 
sirens have luminosity errors ctqw ~ 1% and number 
density new ~ 100/ (47r sr). On the other hand, GW 
standard sirens can determine 1-point functions of the 
matter density better than other methods, in particu- 
lar the probability distribution of lensing magnification. 
This could be useful for distinguishing between different 
dark matter models 1351. 



V. DISCUSSION 

We have shown that observations of the gravitational 
waves emitted by binary compact object inspirals can be 
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